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We propose a method to quantify the Standard Model uncertainty in B — » Kit decays using the 
experimental data, assuming that power counting provides a reasonable estimate of the subleading 
terms in the l/m& expansion. Using this method, we show that present B — + Kit data are compatible 
with the Standard Model. We analyze the pattern of subleading terms required to reproduce the 
B — > Ktv data and argue that anomalously large subleading terms are not needed. Finally, we find 
that S Ks ^o is fairly insensitive to hadronic uncertainties and obtain the Standard Model estimate 
S K = 0.74 ± 0.04. 



A decade of physics studies at the B factories produced 
the impressive set of results on B —* Kn decays summa- 
rized in Table HI As data became more and more accu- 
rate, phenomenological analyses based on flavour symme- 
tries and/or hadronic models were not able to fully repro- 
duce the data. This led several authors to introduce the 
Kit puzzle in its different incarnations [3, 0] . In partic- 
ular, the difference AA C p = A C p(K + Tr°) - A C p(A + 7r~) 
has recently received considerable attention, following 
the new measurement AAqp = 0.164 ± 0.037 published 
by the Belle collaboration [3j. It has been argued that 
AAcp could be a hint of New Physics (NP), but alterna- 
tive explanations within the Standard Model (SM) have 
also been considered. 

To understand whether B — ► A"7r decays are really puz- 
zling, possibly calling for NP, one has to control the SM 
expectations for the B — * Kir amplitudes with a level of 
accuracy dictated by the size of the potential NP con- 
tributions. Thanks to the progress of theory in the last 



Decay Mode 


HFAG average 


global fit 


fit prediction 


10 6 BR(if+7r _ ) 


19.4 ±0.6 


19.5 ±0.5 


19.7 ± 1.0 


10 6 BR(_ft: + 7r ) 


12.9 ±0.6 


12.7 ±0.5 


12.4 ± 0.7 


10 6 BR(lf 7r + ) 


23.1 ± 1.0 


23.8 ±0.8 


24.9 ± 1.2 


10 6 BR(K { \°) 


9.8 ±0.6 


9.3 ±0.4 


8.7 ±0.6 


A CP (K+7v') [%] 


-9.8 ±1.2 


-9.5 ± 1.2 


3.9 ±6.8 


A CP (K + 7v°) [%} 


5.0 ±2.5 


3.6 ±2.4 


-6.2 ±6.0 


A C p(K°it+) {%] 


0.9 ±2.5 


1.8 ±2.1 


6.2 ±4.5 


C(K s ir°) 


0.01 ±0.10 


0.09 ±0.03 


0.10 ±0.03 


S(K S 7T ) 


0.57 ±0.17 


0.73 ± 0.04 


0.74 ± 0.04 


AA C p [%] 


14.8 ± 2.8 


13.1 ± 2.6 


1.7 ±6.1 



TABLE I: Experimental inputs and fit results for B — > Kn. 
For each observable, we report experimental results (BR cxp 
and A^p) 0, 13, H taken from HFAG @], the results of the 
fit using all the constraints (third column) and the prediction 
obtained using all constraints except the considered observ- 
able (fourth column). For AAcp, the prediction is obtained 
by removing both Acp(K + ir°) and Acp(K + tt~) from the fit. 



few years, we know that two-body non-leptonic B decay 
amplitudes are factorizable in the infinite &-quark mass 
limit, i.e. computable in terms of a reduced set of univer- 
sal non-perturbative parameters 0, H, Q. However, the 
accuracy of the predictions obtained with factorization is 
limited by the uncertainties on the non-perturbative pa- 
rameters on the one hand and by the uncalculable sub- 
leading terms in the l/mj expansion on the other. The 
latter problem is particulary severe for B — ► A7r decays 
where some power-suppressed terms are doubly Cabibbo- 
enhanced with respect to factorizable terms (lo| . In- 
deed factorization typically predicts too small B — > Kir 
branching ratios, albeit with large uncertainties. The in- 
troduction of subleading terms, certainly present at the 
physical value of the b quark mass, produces large ef- 
fects in branching ratios and CP asymmetries, leading to 
a substantial model dependence of the SM predictions. 
Given this situation, NP contributions to B — > A7T am- 
plitudes could be easily misidentified. 

In this paper, we suggest a method to estimate the SM 
uncertainty given the experimental data, assuming that 
subleading terms are at most of order 1 / m b . 1 This pro- 
cedure provides a solid starting point for NP searches. 
Clearly, we are not sensitive to the presence of NP con- 
tributions of the same size as the subleading corrections 
to factorization. 

We now describe our method in detail. We start with 
a general parametrization of the B — > A'7r amplitudes 
derived from the one in ref. The decay amplitudes 
are given by: 

A(B+^K°ir + ) = -V ts V t * b P + V us V: b A, 
A(B+^K+7r°) = l=(V u V* b {P + APx + AP 2 ) - 

V us V: b (E 1 +E 2 + A)), 
A(B°^K+tt-) = V ts V* b {P + APi) - VusV^Ex , 



1 An early attempt at this method was presented in ref. 111! . 
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FIG. 1: P.d.f. obtained from the global fit for AA C p (left) 
and for S(K s tt°) (right). 



- — (V ts V t l(P - AP 2 ) 

v us v: b E 2 ) . 



(1) 



In terms of the parameters of ref. [121 ] . our parameters 
read 





= Ex(s,q,q;B,K,%) 


-Pf lM (s,q;B,K,ir), 


E-2 


= E 2 (q,q,s;B,ir,K) 


+ P? IM (s,q;B,K,ir), 


A 


= A 1 (s,q,q;B,K,ir) 


-P? IM ( S ,q;B,K,Tr), 


P 


= Pi(s,d;B,K,n), 




APi 


= P x (s,u;B,K,-k)- 


Px{s,d;B,K,Tr), 


AP 2 


= P 2 (s,u;B,ir,K)- 


P 2 (s,d;B,TT,K). 



(2) 



With respect to the most general parametrization, wc 
have neglected isospin breaking in the hadronic matrix el- 
ements of the effective weak Hamiltonian, yet fully retain- 
ing the effects of the electroweak penguins (EWP). This 
assumption reduces the number of independent param- 
eters and removes the dependence on meson charges in 
the arguments of the parameters on the r.h.s. of eqs. ©, 
where q denotes the light quarks. 




Our procedure is to fit the hadronic parameters to the 
experimental data, taking into account the hierarchy be- 
tween leading and subleading terms in the 1/m^ expan- 
sion by imposing an upper bound to subleading correc- 
tions. Only the correction to the dominant penguin am- 
plitude is well determined by the fit. The information 
on the subdominant terms is limited, while their pres- 
ence contributes to the theoretical uncertainty. The the- 
oretical error on the predicted observables is thus deter- 
mined by the allowed range for the subleading parame- 
ters. While quantifying this range is somewhat arbitrary, 
extreme situations in which the leading and subleading 
terms are comparable would imply a failure of the infi- 
nite mass limit. Of course, one has to be careful about 
possible parametric or dynamical enhancements which 
could invalidate the power counting. Chirally-enhanced 
terms in B — ► Kir amplitudes are well-known examples 
of terms that are formally subleading but numerically 
of 0(1). We have therefore included them in the lead- 
ing factorized amplitudes. We now quantify the allowed 
ranges we use for subleading corrections. To this aim, we 
write each parameter as follows: 



Ey 


= e! + 


Fr(Ei) , 


E 2 


= eU 


F r(E 2 )e lS{E2) , 


A 


= A F + 


Fr{A)e lS{A) , 


P 


= P F + 


Fr{P)e lS{p) , 


APr 


= APf 


+ Pa cm r(AP 1 )e j5(APl) , 


AP 2 


= AP 2 F 


+ Fa cm r(AP 2 )e l5 ^\ 



(3) 



where the factorized amplitudes in the limit — > oo are 



= A nK — a.\ — eft + a\ 



H,EW 



t 4,EW I > 



Aktt [ -012 - g ( a 3,EW ~ a l.Ew\ 



+AnK ( 0% - a l - - a A,EW, 



A F 

pF 

APf = 
AP F = 



A nK ( -a" + "4 + 2 ( a 4,EW ~ ' ! 



4,EW. 



A-xK ( —"4 + 7^ a \,EW j ■ 

— AirK 2 a 4-EW ' 

, 3 c 
Ak-k 2 a 3, EW ' 



(4) 



in terms of the parameters a defined in eq. (31) of 
ref. [l3| • We note that we have discarded non-factorizable 
contributions to the chirally enhanced terms. Further- 
more, 



FIG. 2: P.d.f. obtained from the global fit for the parameter 
r(Ei) defined in eqs. ©. 



A n K = G F /V2m|/^P T (0), 
A k , = G F /V2m 2 B UF K {0). 



(5) 
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FIG. 3: ID and 2D p.d.f.'s obtained from the global fit for the parameters r(E2), 8{E2), r(P), S(P), and r(A), 3(A) defined in 
eqs. ©. 



The coefficient F in eqs. © sets the normalization of 
subleading corrections and is equal to A^k computed 
using the central value of the form factor. The phase 
convention is chosen such that the power correction to 
E\ is real. 

The subleading terms in units of F are given by 
r{X) = [0,0.5] for X = {E x , E 2 , A, APi, AP 2 }. Since 
r{P) is very well determined by the fit, for computa- 
tional efficiency we used r{P) = [0,0.2]. For the sake of 
comparison, Ref. [l]| quotes a value of 0.09l~QQg for the 
contribution to r(P) from penguin annihilation, compat- 



ible with the range we use. All strong phases vary in the 
range [— 7r, it]. 

Using the ranges above for the hadronic parameters 
and the input parameters reported in Table |TTl we per- 
form a fit to the data in Table U using the method de- 
scribed in ref. Flat priors are used for the hadronic 
parameters. Two sets of results are summarized in Ta- 
ble HI On one hand, when using all the experimental 
information as input we test the consistency of the SM 
description of the decay amplitudes in a global fit. On the 
other hand, by removing one of the inputs from the fit 
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FIG. 4: ID and 2D p.d.f.'s obtained from the global fit for the parameters r(APi), 6(APi) and r(AP 2 ), S(AP 2 ) defined in 
eqs. ©. 





0.1307 GeV 
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0.1598 GeV 




0.27 ±0.08 


pB^K jpB^n 


1.20 ±0.10 




1.546 ■ 10~ 12 ps 


T~B+ 


1.674 • 10~ 12 ps 


m,B 


5.2794 GeV 


fB 


0.189 ±0.027 GeV 


m-r: 


0.14 GeV 


rriK 


0.493677 GeV 


A 


0.2258 ±0.0014 


A 


0.810 ±0.011 


p 


0.154 ±0.022 


f) 


0.342 ±0.014 



TABLE II: Input values used in the analysis. Form factors are 
taken from lattice QCD calculations }14| . CKM parameters 
have been taken from ref. [Iff. Wave function parameters can 
be found in Table 1 of ref. flal. 



we obtain a prediction of the corresponding experimen- 
tal observable, using all the other inputs to constrain the 
hadronic parameters. 

Two main results are obtained from the global fit: i) 
the BR values are well reproduced, and they are fairly 
insensitive to the l/mj contributions, but for the CKM- 
enhanced charming penguin P. ii) The values of the Aqp 
are well reproduced, thanks to the l/m& contributions. 
In particular, the presence of AP2 (E2 +A) in the CKM- 
enhanced (CKM-suppressed) part of the B + — > K + n a 



amplitude (see Eq. (TT])) allows to obtain simultaneously 
a positive value of Acp(K + tt°) and a negative value of 
A C p(K + Tr~). This is shown in the left plot of Fig. [U 
where the output distribution of AAqp is fully consistent 
with the experimental world average AAqp = 0.148 ± 
0.028. 

The results for the hadronic parameters are shown in 
Figs. OB] Both the charming penguin parameters r(P) 
and 5(P) are well determined, in agreement with the old 
results of ref. In particular, r(P) is found to be of 
0(1/ mb), as expected from the power expansion in QCD 
factorization. Small values of r(A) are favoured, although 
values as large as 0.5 are not excluded. However, a large 
r(A) requires a 5(A) small and negative. The correc- 
tions to the parameter E2, on the other hand, are pushed 
towards the upper half of the allowed range, namely 
0.3 ± 0.5, showing a preference for a large correction to 
the color-suppressed emission amplitude [1, 03, HH • How- 
ever, we have checked that the p.d.f. for r(E2) falls for 
values larger than 0.6 (although there are other allowed 
regions for r(E2) S> 1, see below). No information on the 
other parameters can be extracted from the fit, but for a 
slight modulation of the phases in the region of absolute 
values close to the upper bound. 
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FIG. 5: P.d.f. obtained from the global fit for Im(r) vs. Re(r) 



We have checked that the result for (APi + AP 2 ) / (Pi + 
E2) is in agreement with the prediction of ref. [l9( (ob- 
tained in the SU(3) limit neglecting left-right electroweak 
penguins). To quantify this statement, we define, follow- 
ing ref. [20j, the SU(3) breaking ratio of matrix elements 



(K*{I = SI2)\Q-\B) 
(Ktt(I = 3/2)\Q+\B) 



(6) 



In factorization, this ratio is tiny due to the fact that 



fkF 1 



}kF* 



so that r - f K F B ^^+f 7r F t 
O(10~ 2 ). However, this cancellation is not related to 
SU(3) (in fact, it also holds for B — > K*ir, where the 
SU(3) argument does not apply). More generally, one 
expects \r\ < 20%. In Fig. [5] we present the value of r 
obtained from our global fit, yielding \r\ = 0.20 ± 0.08. 
The fit is fully compatible with the general expectations 
on SU(3) breaking. The factorization predictions are also 
compatible with the fit result, although the fit prefers 
larger values of SU(3) breaking. 

Going back to the parameters on the r.h.s. of Eq. ([3]), 
we can conclude that pp IM is not the dominant source 
of power corrections in E\ , E2 and A as this would imply 
definite correlations among E\, E% and A which are not 
observed. 

Another mechanism for reproducing the Kir data pro- 



posed in the literature HI, 21, 2^] is a NP contribution 



enhancing the EWP amplitudes with a new CP-violating 
phase. While we do not include NP phases in our anal- 
ysis, we checked that removing subleading corrections to 
emissions and annihilations and allowing r(AP 2 ) to vio- 
late the 1 / raj, power counting, it is not possible to repro- 
duce the Kir data. 

The predictions for the BR, obtained by removing them 
one by one from the fit, show that the observed values 
can be easily explained, all the values being in the ±1ct 
range, the error on the prediction being comparable to 
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FIG. 6: Compatibility plots for ^4cp(^' + 7r ) (upper) and 
Acp(K + tt~) (lower). The cross denotes the experimental val- 
ues. The colour code indicates the level of compatibility with 
the SM prediction. 



the experimental one. On the other hand, the error on 
the predictions for Aqp is much larger than the exper- 
imental precision (up to a factor six for Aqp(K + tt~)). 
Within these large uncertainties, the predictions are in 
agreement with the experimental values at the l-2a level, 
as shown in Fig. [S] 

The choice of the upper limit for the subleading terms 
used in our fit clearly dictates the theoretical error as- 
sociated to the fit predictions. For example, raising the 
upper limit from 0.5 to 1 the error on the fit prediction 
for AAcp increases from 0.06 to 0.09. On the other hand, 
the results of the global fit are fairly independent of this 
choice provided that the upper limit is large enough, as 
shown in Fig. [71 In fact, our point is that a good fit 
of the experimental data can be obtained for subleading 
terms compatible with power counting. Once a good fit 
is obtained, the dependence on the upper bound becomes 
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FIG. 7: Some fit results as functions of the upper bound on power corrections. 
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out using the experimental information on both observables. 
The cross represents the experimental result. 



negligible. On the other hand too small values of the up- 
per limit would result in a worse agreement between the 
theory and the data, showing that the factorization for- 
mulae need to be completed with non-perturbative 1 / m& 
corrections to give a good description of the data. 

Removing both S^ v0 and Cj^^o from the fit, an in- 
teresting prediction can still be obtained for the param- 
eters of the B° — > Kstt time-dependent CP asymmetry. 
We get Ck s7T ° = 0.10 ± 0.04, in good agreement with the 
experimental measurement, and Sk s ti° = 0.74 ± 0.04, 
which is compatible with the experimental world aver- 
age at the ~ la level. In Fig. [8] we show the selected 
region on the Sx s7r o-CK S TT° plane, compared to the ex- 
perimental determination. Both the prediction and the 
measurement are limited by the experimental precision, 
since the other Kir data are a crucial ingredient in our fit. 
For example, reducing all experimental errors by a factor 



of two, the error on the fit prediction for S Ks7T o decreases 
to 0.03, while the error on C Ks ^o decreases to 0.02. It 
is then mandatory to improve the experimental informa- 
tion. Considering the difficulties related to the study of 
B° — > K°tt° in the crowded environment of LHC, Su- 
perB [23| appears as the ideal facility to accomplish this 
task. 



Recently, ref. 
and C 



K s tt 



24j pointed out a correlation between 
Using the experimental value of 



+0.01+0.00+0.00 Q- 

-0.07-0.10-0.06- 01 



C Ks7r o they obtained S Ks7r o = 0.99 
ilar results were found in ref. (25|. Both papers make 
the following assumptions: AI = 3/2 amplitude fixed 
from 7T7T data using SU (3) symmetry (neglecting also left- 
right electroweak penguins). Under these assumptions, 
they solve for the amplitudes A 00 = A(B° — > _ftf 7r°), 

A_| = A(B° — ► K + ir~) and the CP-conjugate ones 

Aqq^ , up to a four- fold ambiguity. This ambiguity 

can be lifted using phenomenological arguments partly 
based on SU (3) and involving charged B — > Ktt modes, 
further neglecting annihilations. Both papers find a 
large value of 0oo = ar g(A)o^oo) ~ 42° leading to a 
prediction for Sk sw o close to one. We have repeated 
the anal ysis and were able to reproduce the results of 
refs. 24, 25|. In addition, we computed the values of the 



relevant hadronic parameters (in our notation: and 
P) corresponding to the four solutions for the amplitudes 
Aoo,A.f ,Aqq,A-i In particular, neglecting annihila- 
tions and AI = 1/2 EWP, the solution with O o ~ 42° 
has a value of P giving a BR(B+ -> K °7r+) - 18 x 10~ 6 , 
incompatible with the measured value (23. 1± 1.0) x 10~ 6 . 
In any case, one gets a huge value of \E2/E\\: using the 
input of ref. [H, we find E 2 /E 1 = 1.9e- a76 ° . Clearly 
this value is not compatible with factorization and would 
imply a breakdown of the heavy quark expansion. In our 
fit, by limiting the range of the power corrections, wc 
discarded this possibility. In fact, we have shown that 
a good agreement with the experimental data is possi- 
ble without introducing huge corrections to factorization. 
Another recent analysis, presented in ref. 26], obtained 
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a good agreement with experimental data, fixing the ra- 
tio of EWP to current-current operator matrix elements 
using QCD factorization and fitting all other matrix ele- 
ments. The range found in ref. [H for \E 2 /E t \ = [0.52,3] 
can possibly be compatible with both our findings and 



the results of refs. [2J, |25| . Indeed it may overlap with 
our findings in the lower r ang e allowed for \Ei/Ei\ but 
also with those of refs. [Hldll in the upper range where 
lE^/Ei] violates the power counting. 

In this letter, we presented a data-driven method 
to estimate the hadronic uncertainties in B — > Kir 
amplitudes compatible with the l/m& expansion. This is 
a basic requirement to meaningfully look for NP in these 
channels. We found that Kir data can be accounted for 
by the SM, including direct CP violation. CP violating 
asymmetries are predicted with a large uncertainty, 
except for Sx s7r o and C^ s7r o, where the theoretical 
error is much smaller than the experimental one. Thus, 
these asymmetries are a better place to look for NP 
than direct CP violation in the other B — > Kit decay 
modes, where possible NP contributions are obscured by 
hadronic uncertainties. 

Note added 

During the completion of this work, we were informed 
that similar results have been obtained by M. Duraisamy 
and A. Kagan in an ongoing analysis of power corrections 
to B — > PP, PV, and VV decays. Earlier results by the 
same group can be found in ref. [27l |. 

We thank R. Fleischer for useful discussions. We ac- 
knowledge partial support from RTN European contracts 
MRTN-CT-2006-035482 "FLAVIAnet" and MRTN-CT- 
2006-035505 "Heptools". M.C. is associated to the Di- 
partimento di Fisica, Universita di Roma Tre. E.F. and 
L.S. are associated to the Dipartimento di Fisica, Uni- 
versita di Roma "La Sapienza" . 



[1] 



J. Buras and R. Fleischer, Eur. Phys. J. C 16 (2000) 
M. Gronau and J. L. Ros- 



[2] 

[3] 
[4] 



97 arXiv:hep-ph/0003323 
ner, Phys. Lett. B 572 (2003) 43 [arXiv:hep-ph/0307095| ; 
A. J. Buras, R. Fleischer, S. Recksiegel and F. Schwab, 
Eur. Phys. J. C 32 (200 3) 45 [arXiv:hep-ph/0309012| ; 
R. Fleischer, |arXiv:hep-ph/0701217| A. Jain, I . Z. Roth- 
stein and I. W. Stewart, arXiv:0706.3399 [hep-ph]; 
S. Baek, Phys. Lett. B 659 (2008) 265 arXi v:0707.2838 
[hep-ph]]; C. S. Ki m, S. Oh and Y. W . Yoon, Phys. Lett. 
B 665 (2008) 231 |arXiv:0707.2967l [hep-ph]]. 
S. Baek and D. London , Phys. Lett. B 653 (2007) 249 
[arXiv:hep-ph/0701181] . 

[The Belle Collaboration], Nature 452 (2008) 332. 
K. Abe et al. [Belle Collaboration], Phys. Rev. D 76 
(2007) 091103 [arXiv:hep-ex/0609006j ; B. Aubert et al. 
[BABAR C ollabo ration], Phys. Rev. D 77 (2008) 012003 
arXi v:0707.2980l [hep-ex]] . 



[5] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 
75 (2007) 012008 [arXiv:hep-ex/060 8003 ; B. Aubert et 
al. [BABA R Collaboration], Ph ys. RevTLett. 97 (2006) 
171805 [arXivihep-ex /0608036] ; K. Abe et al, Phys . 
Rev. Lett. 99 (2007) 121601 [arXiv:hep-ex/0609015] ; 
M. Morello [CDF Collaboration] , Nucl. P hys. P"roc. 
Suppl. 170 (2007) 39 [arXiv:hep-ex/0612018| ; B. Aubert 
et al. [BABAR Collaboration], Phy s. "Rev. Lett. 99 
(2007) 021603 [arXiv:hep-ex/0703016] ; B. Aubert et al. 
[BABAR Collabo ration], Phys. Rev. D 76 (2007) 091102 
|arXiv:0707.2798l [ hep-ex]]. 
E. 



[6] 
[7] 



[8] 



[9] 



[10] 



[11] 
[12] 
[13] 
[14] 

[15] 
[16] 



Barberio et al., arXiv:0808.1297 and online update at 
http://www.slac.stanford.edu/xorg/hfag 
M Beneke, U. Buchalla, M Neubert and 
C. T. Sachrajda, Phy s. Rev. Lett. 83 (1999) 1914 
[arXiv:hep ph/9905312j ; M. Beneke, G. Buchalla, 
M. Neubert and C T. Sachrajda, N ucl. Phys. B 
591 (2000) 313 [arXiv:hep-ph/0006124| ; M. Beneke 
and M. Neubert, Nu cl. Phys. B 651 (2003) 225 
[arXiv:hep-ph/021008"5] ; M. Beneke and S. Jager , 
Nucl. Phys. B 751 (2006) 160 [arXiv:hep-ph/0512351| ; 
M. Beneke and S. Jage r, Nucl. Phys. B 768 (2007) 51 
[arXiv:hep-ph/0610322] . 

H. n. Li and H. L. Yu Phys. Rev. Lett. 74 (1995) 4388 
[arXiv:hep-ph/9409313j ; H. N. Li and H. L. Yu, Phys. 
Lett. B 353 (1995) 30 1; H. n. Li and H. L. Y u, Phys. 
Rev. D 53 (1996) 2480 [arXiv:hep-ph/9411308j . 
C. W. Bauer, D. Pirjol and I. W. Stewart , Phys. Rev. D 
65 (2002) 054022 [arXiv:hep-ph/0109045] ; C. W. Bauer, 
S. Fleming, D. Pirjol, I. Z. Rothstein and I. W. Stewart, 
Phys. Rev. D 66 (2002) 014017 [arXiv:hep-ph/0202088] ; 
C. W. Bauer, D. Pirjol and I. W. Stewart, Phys . 
Rev. Lett. 87 (2001) 201806 [arXiv: hep-ph /0 107002] ; 
J. g. Chay and C Kim Phys. Rev. D 68 (2003) 071502 
[arXiv:hep-ph /0301055 j ; J. Chay and C Kim , Nucl. 
Phys. B 680 (2004) 302 [arXiv:hep-ph/0301262| . 
M. Ciuchini, E. Franco, G. Martinelli and L. Silvestrini, 
Nucl. Phys. B 501 (1997) 271 [arXiv:hep-ph/9703353] ; 
M. Ciuchini, R. Contino, E. Franco, G. Martinelli and 
L. Silvestrini, Nucl. P hys. B 512 (19 98) 3 [ Erratum- 
ibid. B 531 (1998) 656] [arXiv:hep-ph/9708222| ; M. Ciu- 
chini, R. Contino, E. Franco, G. Martinelli and 
L. Silvestrini, Nucl. In strum. Meth. A 408 (1998) 28 
[arXiv:hep-ph/9801420"]; P. Zenczykowski, Phys. Lett. 
B 590 (2004) 63 [arXiv:hep-ph/0402290,i ; M. Ciu- 
chini, E. Franco, G. Martinelli, A. Masiero, M. Pierini 
and L. Silvestrini, |arXiv:hep-ph/0407073"| M. Ciuchini, 
E. Franco, G. Martinelli, M. Pierini and L. Silvestrini, 
Prepared for 32nd International Conference on High- 
Energy Physics (ICHEP 04), Beijing, China, 16-22 Aug 
2001 

M. Pierini, talk given at CKM 2006, December 12-16, 
2006, Nagoya (Japan), http://tinyurl.com/6kkaxb. 
A J. B uras and L. Silvest rini, Nucl. Phys. B 569 (2000) 
3 [arXiv:hep-ph/9812392] . 

M. Beneke and M. Neubert , Nucl. Phys. B 675 (2003) 
333 [arXiv:hep-ph/0308039l . 

A. Abada, D. Becirevic, P. Boucaud, J. P. Leroy, V. Lu- 
bicz and F. Mescia, N ucl. Phys. B 619 (2001) 565 
[arXiv:hep-lat /001 1065] . 
M. Bona et al. [UTfitCollab oration 
081 [arXiv:hep-ph/0606167l . 
M. Ciuchini et al, JHEP 0107 



arXiv:hep-ph/0012308 



JHEP 0610 (2006) 
(2001) 013 



8 



[17] 



C. W. Chiang, M. Gronau, J. L. Rosner and 

D. A. Suprun, Phys . Rev. D 70 (2004) 034020 
[arXiv:hep-ph/0404073j ; Y. Y. Charng and H. n. Li, 



Phys. Rev. D 71 (2005) 014036 arXiv:hep-ph/0410005 
C. S. Kim, S. Oh and C. Yu, Phy s. Rev. D 72 
(2005) 074005 [arXiv:hep-ph/0505060| ; M. Gronau 
and J. L. Rosner, P hys. Lett. B 644 (2007) 237 
[arXiv:hep-ph/0610227] . [22] 

[18] V. Barger, C. W. Chiang, P. Langack er and H. S. L ee, 
Phys. Lett. B 598 (2004) 218 [arXiv:hep-ph/0406126] . 

[19] M. Neubert and J. L. Rosner, P hys. Lett. B 441 
(1998) 403 [arXiv:hep-ph/9808493] ; M. Neubert and 
J. L. Rosner, Phys. Rev. Lett. 81 (1998) 5076 [23] 
[arXiv:hep-ph79 809311 ; M. Gronau, D. Pirjol and [24] 
T. M. Yan, Phys. Rev. D 60 (1999) 0340 21 [Err atum-ibid. 
D 69 (2004) 119901] [arXiv:hep-ph/9810482j; M. N eu- [25] 
bert, JHEP 9902 (1999) 014 |arXiv:hep-p h/9812396 . 

[20] M. Ciuchini, M. P ierini and L. Silvestrini , Phys. Rev. D [26] 
74 (2006) 051301 [arXiv:hep-ph/0601233| . 

[21] Y. Grossman, M Neubert and A. LT Kagan, JHEP [27] 
9910 (1999) 029 [arXiv:hep-ph/9909297|; T. Yoshikawa , 
Phys. Rev. D 68 (2003) 054023 [arXiv:hep-ph/0306147| ; 
S. Mishima and T. Yoshikawa, Phys. Rev. D 70 (2004) 



094024 arXiv:hep-ph/0408090 ; S. Baek, P. Hamel, 
D. London, A. Datta and D. A. Suprun, Phys. Rev. D 71 
(2005) 057502 [arXiv:hep-ph/0412086| ; R. L. Arnowitt, 
B. D utta, B. Hu and S. Oh , Phys. Lett. B 633 (2006) 
748 [arXiv:hep ph /05092331 ; M. Imbeault, S. Baek 
and D. London , Phys. Lett. B 663 (2008) 410 
[arXiv:0802. 11751 [hep-ph]] . 

A. J. Buras, R. Fleischer, S. Recksiegel and 

F. Schwab, Phys. R ev. Lett. 92 (2004) 101804 

[arXiv:hep-ph/0312259] ; A. J. Buras, R. Fleischer, 

S. R ecksiegel and F. Schwa b, Nucl. Phys. B 697 (2004) 

133 [arXiv:hep-ph/0402112j. 

M. Bona et ffl/., larXiv:0709.0451l [hep-ex]. 

R. Fleischer, S. Jager, D. Pirjol and J. Zupan, Phys. Rev. 

D 78 (2008) 111501 [arXiv: 0806. 29001 [hep-ph]] . 

M. Gronau and J. L. Rosner, Phys. Lett. B 666 (2008) 

467 |arXiv:0807.3080l [hep-ph]]. 

T. Feldma nn, M. Jung and T. Mannel, JHEP 0808 
(2008) 066 [arXiv:0803.3"729l [hep-ph]]. 
A. Kagan, talk given at CKM 2006, December 12-16, 
2006, Nagoya (Japan), http://tinyurl.com/679qkv. 



